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Abstract 

We establish existence of global-in-time weak solutions to the one dimensional, 
compressible Navier-Stokes system for a viscous and heat conducting ideal poly- 
tropic gas (pressure p — K6/t, internal energy e = c v 9), when the viscosity /i is 
constant and the heat conductivity k depends on the temperature 6 according 
to k(8) — RO 13 , with < [3 < §. This choice of degenerate transport coefficients 
is motivated by the kinetic theory of gasses. 

Approximate solutions are generated by a semi-discrete finite element scheme. 
We first formulate sufficient conditions that guarantee convergence to a weak 
solution. The convergence proof relies on weak compactness and convexity, and 
it applies to the more general constitutive relations /i(#) = fi9 a , k(9) = R9 13 , 
with a > 0, < (3 < 2 (a, R constants). We then verify the sufficient conditions 
in the case a = and < [3 < §. The data are assumed to be without vac- 
uum, mass concentrations, or vanishing temperatures, and the same holds for 
the weak solutions. 
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1 Introduction 



In one space dimension (ID) the Navier-Stokes system for compressible flow 
of a heat conducting and viscous fluid takes the following form in Lagrangian 
coordinates: 



n = u x 



£t + (up)* 



(mass conservation) (f.l) 
(momentum balance) (1-2) 

(energy conservation) (f.3) 
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Here :r=Lagrangian space variable, t=time, and the primary dependent variables 
are specific volume r, fluid velocity u, and temperature 9. The specific total 
energy £ = e + \u 2 , where e is the specific internal energy. The pressure p, 
the internal energy e, and the transport coefficients {i (viscosity) and k (heat 
conductivity), are prescribed through constitutive relations as functions of r 
and 9. The thermodynamic variables are related through Gibbs' equation de = 
OdS — pdr, where 5'=specific entropy. See |22| for a derivation of the model. 

In this article we establish existence of weak, global-in-time solutions to the 
field equations {^-(T^-|[L3} for given initial and boundary data. The flow 
domain is a finite interval SI := (0, L) in Lagrangian coordinates. The data we 



consider are standard (see Theorem 2.1 for precise regularity assumptions): r, 



u, 9 are prescribed initially on SI, while it and 9 satisfy homogeneous Dirichlet 
and Neumann conditions, respectively, on dSl; 

u = 0, 9 X =0, onfflx(0,T). (1.4) 

Our main interest concerns the transport coefficients /i and k. Their depen- 
dence on r and 9 will obviously influence the solutions of the field equations as 
well as the mathematical analysis. Even for ID flows there is a wide gap between 
the models furnished by physical theories, and the models covered by a satis- 
factory existence theory. We focus on the case of gases for which kinetic theory 
provides constitutive relations and we consider only ideal, polytropic gases: 

p = p{9,r) = K-, e = c v 9, (1.5) 

T 

where K (specific gas constant) and c v (specific heat at constant volume) are 
positive constants. We scale c v to unity such that £ = 9 + \u 2 . 

According to the first level of approximation in kinetic theory the viscosity 
fi and heat conductivity k are functions of temperature alone. Furthermore, 
the functional dependence is the same for both coefficients. See Chapman & 
Cowling [5] or Vincenti & Kruger 24 for a thorough discussion of these issues. 



If the intermolecular potential varies as r a , r=intermolecular distance, then /1 
and k are both proportional to a certain power of the temperature: 

/i, K (X 9 2a . 

For Maxwellian molecules (a = 4) the dependence is linear, while for elastic 
spheres (a — > +00) the dependence is like \[9. In any case 

/i = fi9 b and k = R9 b for some b £ (|,+oo), (1-6) 

where fl and R are constants. In particular, the transport coefficients tend to zero 
with 9. The discrepancy mentioned above is illustrated by the fact that, beyond 



the regime of small and sufficiently smooth data 13 , there is no global-in-time 
existence result currently available for the Navier-Stokes model ( |l.l[ )- |L2l )-( 1.3 1, 



with constitutive relations (1.51 and (1.6 1. 

To put our main result Theorem |2.1| in perspective let us contrast the last 
statement with what is known about existence of solutions to the compressible 
Navier-Stokes system. (For the purpose of this introduction we concentrate on 
ID flows and we consider only a small selection of the very extensive literature.) 
The seminal work of Kazhikhov & Shelukhin 16 treats the full one-dimensional 
Navier-Stokes system ( 1 . 1 1- ( 1.2 1- ( 1 .3 1 with \\.b\ and constant transport coeffi- 



cients. Building on earlier work by Nash [21], Kanel [12 , and Kazhikhov 15 



global existence and uniqueness of smooth (i.e. W ' ) solutions are established 



in 16 for arbitrarily large and smooth data. A key ingredient in the proof is 
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the pointwise a priori estimates on the specific volume which guarantee that no 
vacuum nor concentration of mass occur. 

Much effort has been invested in generalizing this approach to other cases, 



and in particular to models satisfying (1.61. This has proved to be challenging 



Temperature dependence of the viscosity /i has turned out to be especially prob- 
lematic. On the other hand, one has been able to incorporate various forms of 
density dependence in and also temperature dependence in k. Dafermos [7], 
Dafermos & Hsiao [8] considered certain classes of solid-like materials in which 
the viscosity and/or the heat conductivity depend on density, and where the heat 
conductivity may depend on temperature. However, the latter is assumed to be 



bounded as well as uniformly bounded away from zero. Kawohl 14 considered 
a gas model that incorporates real-gas effects that occur in high-temperature 
regimes. In [l4] the viscosity depends only on density (or is constant) and 
it is uniformly bounded away from zero, while the thermal conductivity may 
depend on both density and temperature. For example, one of the assump- 
tions in [14] is that there are constants kq, Ki > such that k(t,9) satisfies 
Ko(l + 9 q ) < k(t, 9) < Ko(l + 9 q ) where q > 2. This type of temperature depen- 
dence is motivated by experimental results for gases at very high temperatures, 
see Zel'dovich & Raizer [28] . None of these results cover the case of a degenerate 
heat conductivity k(9) = R9 b . 

In the case of isentropic flow a temperature dependence in the viscosity 
translates into a density dependence. For some representative works in this 



direction see 18 , [191, 20 , 25 , [26], 27 , and references therein. Finally, for the 



multi-dimensional Navier-Stokes equations there has recently been established 
various existence results where the viscosity and/or heat conductivity depends 
on t or 9, see [3] and [10| . These models do not cover the case of a gas with 



constitutive relations (II. 51 and (1.61 



Outline The overall approach in the proof of Theorem |2.1| is standard: apri- 
ori pointwise estimates on specific volume and temperature are coupled with 
higher-order integral estimates to provide sufficient compactness to pass to the 
limit in an approximation scheme. However, in both parts of the analysis the 
temperature dependence in the transport coefficients raises some new issues. 

To generate approximate solutions we use a semi-discrete finite element 
scheme. This provides an easy proof of well-posedness of the scheme (Section^ 
while avoiding some of the cumbersome notation of finite difference schemes. In 
Section [4] we formulate certain apriori bounds on the approximations and show 



that these are sufficient for convergence to a weak solution (Theorem 4.1 1. The 
proof employs weak convergence and convexity techniques a la Lions [17] and 



Feireisl 10 . The particular proof we use, working entirely in the Lagrangian 
frame, seems to be new. The argument applies to ideal polytropic gases with 
/i oc 9 a with a > 0, and k oc 9 13 with (3 £ [0,2), thus including the "standard" 
case of constant transport coefficients, as well as a full range of powers pre- 



dicted by kinetic theory (see (1.6 1). More general constitutive relations could 
presumably be included at the expense of more detailed growth conditions. 

In Section [5] we verify the sufficient apriori bounds in the case where k(9) = 
k9® , with /3 6 [0, |), and fi is constant. Again the result covers constant heat 
conductivity and a range of power laws suggested by kinetic theory. However, we 
have not been able to establish sufficient apriori estimates in the case where also 



the viscosity depends on temperature. (We point out in the proof of Lemma 6.1 
where constant viscosity seems to be essential for the argument.) To establish 
the apriori bounds we derive both pointwise estimates and certain energy esti- 
mates. To treat the latter we follow Hoff jTT] and define energy functionals that 
monitor certain weighted _H' 1 -norms in the solution. At this point the tempera- 
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ture dependence in the heat conductivity requires a careful choice of functionals. 
On the other hand, the assumption of constant viscosity allows us to adopt a 
standard argument, with only minor changes, to obtain the necessary pointwise 
estimates on r and 8. For completeness this part is included in Section [6] 



2 Main result 

It is convenient to formulate the approximation scheme for the temperature 
field instead of the total energy £ . Consequently we consider a weak form of the 
temperature equation: 



= M— -pu x . (2.1) 

T I T 

/ x 

In the following definition = (0, L), with < L < oo, and we assume that u~, 
K, and p are given, smooth, non-negative functions of (9,t). 

Definition 2.1 (Weak solution) We say that (r,u,6) is a weak solution of 
the compressible Navier-Stokes system ( |1.1[ ), | |1.2[ ), < |2 . 1 [ ) on Q, x [0,T), with 
boundary conditions (1.4 I and initial data (tq,Uo,0o) satisfying 

(r ,u ,e ) G L 2 {Q) x L 2 {Q) x L 2 {fl) , (2.2) 

provided that r G W 1,2 (0, T; L 2 (Q.)), u € L°°(0, T; L 2 (Q.)) n L 2 (0, T; Wq' 2 (£1)), 
6 e L°° {0,T; L 2 (n)) n L 2 (0,T; W 1 ' 2 ^)), 

Tt = u x , a.e in Q x (0, T), t(x, 0) = tq(x), a.e x G Q, (2-3) 

and that for all <f> G C" c °°(^ x [0,T)): 

u(j>t — \ — - — pi <j>x dxdt + / u 4>odx = Q, (2.4) 
o Jn L r J J n 

f f e<j>t--e x <t>x + \_—-p\u x ct>dxdt+ j e (j) dx = Q. (2.5) 

In addition all terms in (2.4 1 and (2.5 I are required to be mtegrable. 

In the remainder of the paper we restrict ourselves to the case of an ideal poly- 
tropic gas, i.e. (1.5 I holds. Furthermore, for concreteness, and motivated by 
kinetic theory, we only consider the case where fi and k, are proportional to 
(possibly different) powers of 9: 

H(6)=fi6 a , k(60 = k6»' 3 , (2.6) 

where p,, k > 0, and a, j3 > are constants. To state the main result we define 
the following functionals: 



A(t) = sup / \u x \ 2 dx + j / \ut\ 2 dxdt, 
s£(o,t)Jn Jo Jn 



B(t)= sup / L(6)dx+ / / \L'(8 h )x\ 2 dxdt, 
se(o,t)Jn Jo Jn 

V(t)= sup I f -\L'(6) X \ 2 dx+ f f K(6)\9 t \ 2 dxdt , 
se(o,t) 2 Jq t J J n 



where 



L{z) := I' ^ K{ n ) d V d£ . (2.7) 



Our main result is: 
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Theorem 2.1 Consider the one- dimensional, compressible Navier-Stokes sys- 
tem ( |1.2[ ), and ( 2. 1 1 for an ideal, polytropic gas (1.5l. Assume that the 



transport coefficients and re satisfy \2.Q\ with a — and < (3 < 
initial data satisfy 



(t , mo, «o) G L (SI) x W L '\Sl) x VT'"(fi), 



and be such that 



C 1 < to(x) < Co , C 1 < #o(a:), /or a. a. x G O, 



(2.8) 



(2.9) 



where Co > is a constant. 

Then, for any finite time T > there exists a weak solution (t, u, 6) on 
SI x [0, T) o/ \1.2\ , and (2.1 1, «nf/i f/ie boundary conditions (1.4 1 and 

initial data (tq,Uo,Qq). Furthermore, there exists C > 0, depending on the 
parameters K, c v , fi, R, (3, the initial data, and T , such that 



C" 1 < r(x,t) <C, 



?(x,t) > 



a.e m!!x (0,T), 



and 



sup [^(t)+B(t)+D(t)] < C. 

'6(0, T) 



(2.10) 



(2.11) 



TTie weafc solution can be obtained as the pointwise a.e. limit in Q x [0,T) o/ 
solutions to the semi-discrete finite element scheme described in Definition \3.1\ 

The proof is detailed in the following sections and summarized in Section [5. 1| 



Remark 2.1 Concerning the regularity of the initial data we note that the result 
in |11] covers discontinuous data under the assumption of BV regularity of to 
anduQ. The scheme we consider (Definition 3.1) is well-defined for (to, uq,6o) G 



L (fi)l ; the higher regularity in (2.8 1 is required to bound the initial values of 



A, B, V. We also mention that (2.101 may be refined to give time-independent 
density bounds by adopting the techniques in [§]. 



3 Finite element scheme 

We define a semi-discrete finite element scheme approximating the compressible 
Navier-Stokes system. We verify the basic bounds for mass, energy and entropy 
for the scheme, and show that the scheme is well-defined globally in time. 



3.1 Approximation scheme 

Let {Eh}h>o be a family of uniform meshes of SI, where h is the mesh size and 
it is assumed that N = ^ G N. For each such h, the vertices of are Xi — ih, 
i — 0, . . . , N . For each interval £ £ E^, Pit (15) denotes the space of polynomials 
of maximal order k on E. On Eh we define the space of piecewise constants, 

Qh(Sl) = {<!>& L 2 (Sl); (t>\ E G Po(-E), V£ G E h } , 

and the space of continuous piecewise linears, 

V h (Sl) = {v G C(Sl);v\ E G Fi(E), WE G E h } . 

Elements in Qh(Sl) are taken to be continuous from the right, as are derivatives 
of elements in Vh(Sl). In order to incorporate Dirichlet boundary conditions we 
define the space 

Vg(Sl) = {v£V h (Sl);v\an = 0} ■ 
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Figure 1: A solution of the scheme at a fixed time. The specific volume r h and 
temperature 9 h are piecewise constant while the velocity u h is piecewise linear. Note 
that the nodes are numbered from while the elements are numbered from 1. 



For Qh(fi) we define the projection operator 11^ : L P (Q) — * Q^(f2), p G [l,oo] 
by 

/ \n%<j>](x) dx= I <t>(x) dx, VEeE h . (3.1) 

J E J E 

For Vh(f2) we define the projection operator II : W 1 ' 2 (Q.) — > Vj,(fi) by 

[n^](^) = ^), » = o,...,jv, 

such that 

{UU) X =I^{4> X ). (3.2) 
For any q G Q/t(fi) we set [g], := ?|.E i+1 — , where i5i denotes the ith element 

Definition 3.1 (Semi-discrete finite element scheme) Let the constitutive 
relations (1.5 I and (2.6 1 /ioZd and fix h > uiit/i JV := ^ £ N. For initial data 
(to,Uo,Oo) satisfying (2.2 1 we set 

(t \4X) := (n? N.nfciuoi.n^po]), (3.3) 

where Tlh : i p (f2) — > F (fi) is t/ie L 2 -projection. For a given time T > we 
determine functions 

(T h ,u h ,8 h ){t) G Qfc(fi) x v£(0) x \4(fi), t G (0,T), 



sttc/i t/iat 



(3.4) 

\/v h G t£(n) (3.5) 



/ Tf (j) h dx= u^.(/> h dx \/<f> h G Qh(f2) 

in in 

jf ^ = -/ n [^^-P(*V h )]i£ dx 

/ t V dx+'-Y. Gi(r h )lL'(e h )W% 
Jn n i=i 

= / [ M(g y 2 -p(flV>£]^ V/6Q h (!]), (3.6) 

where G* i (r' 1 ) = 2(t h \ Ei + T h \ Bi+1 ) and L{z) is given by &2.71. 

3.2 Basic estimates and well-posedness of the scheme 

To show that the scheme is well-defined we recast ( |3.4[ )-( 3.6 1 as an ODE on the 
(finite dimensional) finite element space Qh(fi) x V^JTl) x Qh(fi). Assuming for 



(i 



now that a solution (r , u ,9 ) exists, we show mass and energy conservation, 
and entropy balance, for the scheme. With (f> h = 1 equation (3.4 1 gives 

r h (x,t) dx= t (x) dx . (3-7) 
Jn 

Next, (3.5 1 with v h — u h and (3.61 with tp h = 1 give conservation of energy: 

£(u h ,e h ) ■- f {\\u h \ 2 + e h ){x,t) dx = £{utei). (3.8) 



Finally, |3j} with ^ h = ^ and {31} with <p h := £ give 



d_ 
dt 



log 9 h + K log r h dx = 



.hah ax h 

i=l 



JV-1 

i ^ G i (r")[L'(fl h )] i 



Gi(r")«(^ t ) Mh 



]? (3-9) 



where, for each i = 1, . . . , N — 1, Oh, 9^ lie between 1^ and 9 h \E i+1 (obtained 
from the mean value theorem applied to L'(9) and |, respectively). Thus the 
discrete entropy 



S(T h ,u h ,9 h ) := \\u h \ 2 + 9 h + Kr h - log(0 h ) - Klog(r h ) 



satisfies 



d 
dt 



S[r\u\9*)dx + j^dx + \Y i 



Lemma 3.1 (Well-posed ness of approximation scheme) Under the same 
assumptions as in Definition 3.1 there exists a unique solution (r h ,u ,9 h ) of the 
approximation scheme ( 3.4 1 -( 3.6 1 . The solution is defined for all times t > 
and there is a constant Co, depending on h and (to,Uo,0o)> but independent of 
time, such that 

Co 1 < r h {x,t), 9 h (x,t) < Co Vx e n, vt > o. 

Proof: We define the finite element space S/^fi) := Qh(fi) x V°(fl) x Qh(£l) 
of (r h , u h , 9 h ) triples and equip it with the standard L 2 (Q,) = L 2 (Q) x L 2 (fi) x 
L 2 (fi) norm. Let S£(Cl) denote the subset S h (Sl) <1 { r\ 9 h > 0}. For each 
z = (21,22,23) £ S^fi) we define F(z) to be the unique element in Sh(fi) 
satisfying 



/ (22)^91 da; + 









21 



(23,2l) 



JV-1 

- £ Gi(zi)[L' (*i)Mfo, 



[{Z2)xq3 — (92)1] da; 

(3.11) 



for all q = (gi, 92,93) £ S^fi), where (•,•) denotes the standard inner product 
on L 2 (fl). Linearity with respect to q in (3.111 implies that the approximation 
scheme |3.4l-|3.6l corresponds to the ODE 



z t = F(z), z(0) = Zq = (to,Uo,0q). 



(3.12) 
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For the given initial data Zq we define the compact subset 



(z 2 ,z 3 )dx< I £(uq,9o) dx, 



K{z{])~{z G S+{Q); f Zl dx< f T & dx, f £(z 2 ,z 3 ) dx < f 
1 Jn Jn Jn Jn 

log(z 3 zf) dx > / log (9o(to) k ) dx\ 

J r> J 



It follows from the estimates (3.7 1-13.10 1 that any solution z(t) of the initial 
value problem (3.121 belongs to K(z ) for all times it is defined. In turn this 
gives a time independent (but /i-dependent) upper bound for ||«(f)||ioo(m, as 



well as for \\z 1 1 (t)\\ 



i,°°(n) 



and 



_1 C*)II 



A direct calculation shows that, 

there 



thanks to these bounds, F is uniformly Lipschitz continuous on K(zq) 
is a constant L (depending on h and zj) such that 



\F(z) - F(y)\\ L 2 in) <L\\z-y 



£ 2 (f2)! 



Standard ODE theory now gives existence and uniqueness of the solution z(t) to 



( 3.12) for all times t > 0. I.e., (j h ,u h ,G h ){t) = (zi(t), z 2 (t), z 3 (t)) is the unique 
solution to the approximation scheme ( 3.4 1- ( 3.6 1 with initial data Zq. □ 



4 Convergence of the approximation scheme 

In this section we formulate certain (strong) apriori estimates on solutions of the 



scheme ( 3.4 1- (3.6 1, and we show that these are sufficient to conclude convergence 
to a weak solution of the system (JTTTj) , ( |1.2[ ), ( |2.1[ |. As noted in the introduction, 



this part of the analysis works for more general powers a and /3 (see ( 2.6 1) than 
those assumed in the statement of Theorem 12. II On the other hand we recall 
that we restrict ourselves to ideal polytropic gases. 

Notation 4.1 Given a set of elements {v h } in a normed space (X, \\ ■ ||), we 
write "v h G(, X" to mean that is bounded independently of h. Weak con- 

vergence is denoted by — 1 and weak-* convergence is denoted — s -. An over bar 



denotes weak L 1 -limit, see e.g. (4.151. Subsequences are not relabeled. A zero 



subscript denotes evaluation at time t — 0. 

For later reference we recall the following general results: 

Lemma 4.1 Let {f h }, {g h } be sequences of functions on a measurable subset 
O of n x (0,T). Assume that f h e b L°°\o) with f h -> / a.e. in O, and 
g h £b L x {0) with g h in ^(O). Then f h g h fg in l}(0). 

Lemma 4.2 (M Theorem 2.11]; Let O C R M , M > I, be bounded and open and 
let g: E — > E be continuous and convex. Let {v n } n>1 be a sequence in L 1 ^) such 

that v n — 1 v in Z, 1 (O), g(v„) S i 1 (0) for each n, and g(v„) — 1 g(v) in L l (0). 
Then g(v) < g(v) a.e. on O and g(v) € L l (0). If, in addition, g is strictly 
convex on an open interval (a,b) C E and g(v) — g(v) a.e. on O, then, passing 
to a subsequence if necessary, v n (x) — * v(x) for a.e. x 6 {y £ O | v(y) G (a, 6)}. 



Lemma 4.3 ( 23 Corollary 4]) Let X C B C Y be Banach spaces with X C 
B compactly. Then, for 1 < p < oo, {v : v G L p (0,T;X),v t G L^O^;^)} 
is compactly embedded in L v (0,T; B) . In the case p = oo, for any r > 1, 
{v : v G L°°(Q,T;X),Vt G L r (0,T;Y)} is compactly embedded in C(0,T;B). 

The main result in this section is: 
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Theorem 4.1 Assume the constitutive relations (1.51 and (2.6 1 withO < (3 < 2. 
Let {(t h ,u h ,9 h )} be a sequence of functions constructed according to Definition 
3.1\for t G [0, T], and assume there exists C > 0, independent of h, such that 



[CI] 



cr 1 < r h {x,t), e h {x,t) 

fi(9 h (x,t)) < C 



N-l fl 

[C2] : sup ||^(t)|| L2(S2) + J2 \ 

*6(0,T) i=1 Jo 



V(i,i)e!!x [0,T) 



dt<C 



[C3] 



u 66 L oo (0,T;L 2 (f2))nL 2 (0,T;W r 1 ' 2 (f7)). 



Then, passing to a subsequence if necessary, (r h ,u h ,6 h ) — > (r,u,8) a.e infix 
(0, T), where (t, u, #) is a uieafc solution to the compressible Navier-Stokes system 
fLTj ) - |2l]) in t/ie sense of Definitions. 



This theorem will be a consequence of Lemmas |4.4| - 14.7| 

Remark 4.1 Clearly, [CI] is satisfied for a — u>/n/e it is equivalent to a 
uniform upper bound on 8 h when a > 0. 



By [CI] — [C3] there exist functions r, u, 6 such that 

u h ^u in L°°(0,T;L 2 (n))nL 2 (0,T;W / "o' 2 (fi)) 
r h ^r in W 1,2 (0,T;L 2 (ft)) 
h ^(9 in L°°(0,T;L 2 (n)), 



(4.1) 
(4.2) 
(4.3) 



if necessary by passing to a subsequence. Here (4.1 1 follows from [C3], (4.2| 
follows from T t h = v% and [C3], and follows from [C2]. 

We will show that the triple (t, u, 6) is a weak solution of the Navier-Stokes 
system. In Section |4.1| we prove that the velocity and temperature converge 
strongly by utilizing Lemma [4. 3| Then a convexity argument is used repeatedly 
in Section [4. 2| to show that the specific volume converges strongly. At this point 
it follows that (r,u,8) is a weak solution of ( | 1 . 1 [ > and (1.2 1. Exploiting this 
and the fact that the pressure is an L 2 function, we find that {u^} converges 
strongly. It then follows that the limits satisfy the temperature equation as well, 
and Theorem I4JJ follows. 



4.1 Strong convergence of u h and h 



Then 



Lemma 4.4 Let {r h ,u h ,6 h } be as in TheoremU.l 

ui 66-L 2 (0,T;W- 1,2 (O)), (4.4) 
and consequently, if necessary passing to a subsequence, 

u h ^u in L 2 (0,T;L 2 (fi)). (4.5) 
Proof: Fix v G i 2 (0,T; Wq' 2 (Q.)) and define functions v h G Vfi(fl) by requiring 



vw h dxdt V«j A 6 Vfc(n) . 



i.e., v h is the L 2 -projection of v onto V£(Q,). Then, by E], there exists a constant 
C > 0, independent of h, such that 



\\V Hi2(0,T;W 1 . a (£l)) < C|kllL 2 (0,T;Wl.2(n))- 
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Using v h as test function in the momentum scheme (3.5 I gives 

rT 









/ Ut v dxdt 




L 


Jo Jn 







UfV h dxdt 



< ll u ' 1 ||L 2 (0,T;W 1 . 2 (n)) 

< C\\ V IU 2 (0,T;lVl. 2 (n))! 



L 2 (0,T;L 2 (fi)) 



+ l|p(^'\T' I )llL 2 (0,T;L 2 (n)) 



where we have used [CI] - [C3]. As v is arbitrary, we conclude (4.4 1. 



In the current situation where both [C3] and (4.4 1 holds, Lemma 4.3 with 
X = W 1,2 (ty, B = L 2 (n), and Y = W~ 1>2 (Q), yields ([Di}. □ 

To analyze the convergence of 8 h we define a dual mesh E^ with vertices 



,N- 



0, xi 



L. 



We define finite element spaces Vf^~(Q), Q^(Q) and corresponding projections 

. 1 

* = nX^[4>]. (4.6) 



II and 11^ exactly as we did for E^. To shorten notation we write 



Then, 



Lemma 4.5 Let {(r h , u h , 6 h )} be as in TheoremU.l 

¥ e b L 2 (0,T; W 1,2 (Q)), 0t £& L 1 (0, T; W _1,1 (f2)) , (4.7) 

and 9 h — *■ in L p (0, T; L V (Q,)), Vp G [1, 4), i/ necessary passing to a subsequence. 
Proof: The dual mesh is defined such that: 



x G (Xi ,x i+1 ), i = 1, 



, N - 



and by [C2] there is a constant C > 0, independent of h, such that 



l h x\\L 2 (Q,T;L 2 (n)) — 



Y l ^dt<c. 



(4.8) 



Since ||0>»(t) - h (i)|| L 2 (n) < C/i||^(t)|| i3(n) , we thus have 

H6I' 1 - 6»' 1 || L 2 (0iT;i 2 (n)) < Cft||6'' 1 :E || L 2 (0iT;I ,2 ( n )) <C7i, 
and [C2] gives 

\\8 h \\L 2 (0,T;L 2 (n)) < ll^ h _ ^ h |ll, 2 (0,T : L 2 (n)) + II ^ II L 2 (0,T;L 2 (fj)) < C(l + h ) : 



which proves (4.7^. For later reference we note that, again by (4.8|, 6 h Gb 
L 2 (0,T;L°°(fi)), whence the same holds for 9 h . Together with [C2] and the 
Cauchy-Schwarz inequality this gives 

6 h G b L 4 (0,T;L 4 (r!)), (4.9) 

which in turn shows that K,(9 h ) Gb L 2 (0, T; L 2 (Q)), since /3 < 2. 

For (4.7l 2 we fix t/j G W 1,00 (f2) and define the discrete effective viscous 

flux by T h := -p(6 h ,T h ). By [CI] - [C3], and the Cauchy-Schwarz 
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inequality we have T h Gb L 2 (0,T; L 2 (tt)). Using H?ip as test function in the 
temperature approximation (3.6 1, we deduce 

i r i r 1 JV ~ 1 

/ 0ty dx\ = \ T h u h x ^ dx~^J2 Gi(T h )lL'(6 h Mn%m 



rr~h h 
\J- U 



a;llL 1 (n) 



+ lhM 



(4.10) 



< IIV-II 



li 2 (n)ll M ^lli' 



JV-l 

!( n) + qk(^)|| L2( n)[E X 1 



where we have used [CI] and that L"(z) — k(z). Integrating (4.10 1 and applying 
[C2]-[C3], together with the bounds on T h and k(6 h ), then gives ET7l) 2 . 

Next we establish compactness of {6 } with respect to the spatial variable. 
Let £ € (— L, L) be arbitrary and define the domain Q.^ = (|£|, L— 1£[). If ft < |£|, 
fix any x £ Q^, let _Ei be the element containing x, and let Ej be the element 
containing x — £. Then 

I e h Or, t) - e h (x ~ e, t) 1 2 = I s' 1 (t) k - 0' 1 (t) I E . 1 2 

Integrating over £7^ x (0,T), and using h < |£| together with [C2], we obtain 

\\e h {x,i)-6 h {x-t,t)\\v( ,Wn t )) < C \^- ( 4 - n ) 
On the other hand, if |£| < ft, then for each < t < T we have (for ft < 1) 



2 


J-1 " 




[v 1 «~ 

^ ft 

.fc=l 






< 


E ft 

.fe=i . 




< C(\£\+h) 


^ ft 

- fc=l 



? h (x,t) -6 h (x-£,t)\ 2 dxdt 



J o i= i 



T»-l 170/1112 



E 



whence [C2] shows that (4.11| holds also for |<J| < ft. From the translation 
estimate (|4.11[), it follows that 6 h € b L 2 (0, T; W2" f,2 (Q)), see m. 



Consequently, Lemma 
Y" = W -1 ' 1 ^), yields 



4.3 



with X = H/2- E > 2 (!T2), e > 0, S = L 2 (ft), and 



6T-S-0 in L 2 (0,T;L 2 (r2)) 



Finally, (4.91 and (4.121 concludes the proof. 
We observe that ( |4.12[ ) and [CI] yields: 



in L p (fix(0,T)) Vl<p<(50. 



(4.12) 

□ 

(4.13) 



4.2 Strong convergence of 



We first note that r t h = a.e. by (3.4 1, whence (4.1 1 and (4.2 1 show that r, u 
satisfy (2.3 1. Passing to the limit a.e. in the other equations requires strong 
convergence of t . We begin by observing that, due to [CI], we may divide by 
r in (2.31 to obtain the " renormalized" equation (used below in (4.191) 



(log r) t = 



a.e. in fi x (0,T) . 



(4.14) 



We proceed to adapt a simplified version of the convexity arguments in 17 , [10] 
to the Lagrangian setting and obtain the required strong convergence. 
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4.1 



Then, 



Lemma 4.6 Let {t\ u h , 6 h ) be as in Theorem 

r h -> t a.e. in Six (0,T) . 
Proof: Step 1: Recalling the notation for weak L 1 -limits we claim that 



pA_«£> a .e. infix (0,T). 
For the proof we fix tp G C™{Sl x (0, T)) n {ip > 0} and set 



(4.15) 



v h {x, t) := J j (y, p| [V' 7 " r T ] (y, i) dy , 

where the integrand is differentiable in time since 

G b L^Til 1 ^)). 
Hence, u£ Gb L 2 (0, T; and by (|47TJ> and ((OJ we have 



(4.16) 







in L\Q,T- L°°(fi)). 



(4.17) 



We then use n^fti ] as test function (3.5 I and integrate in time to obtain 



^dxdt- I I ip[p(6 h ,T h ) -p{6 h ,T)jdxdt 



»{0 h )u h x 



p(9 h 



dx 



o Jn 
1 



+ 



u h Vt dxdt + 



/ / u' 1 life - Vt dxdt —. y^Ji 
Jo Jn 1 J i=1 



(4.18) 



where we have used ( |3.1[) and ( |3.2[ ) . We observe that Ii , J2 — + as h —> 0, due 
to (4.171 and Lemma 4.4 The 73 term satisfies 

\h\ < ||w' 1 || i 2 (0jT . i oc ( n))||n^[u t A ] -vt\\ L2{0T . LHn) y 

where the first factor is bounded by [C3]. A standard interpolation estimate 
[3], together with ( |4.16| l, gives 



|njf [Vt] - Vt I j £2(0,7^1(0)) - ^IK*IU 2 (0,T;Z,1 



(f!)) 



< C7i. 



Recalling (4.131 and Lemma 4.5 and sending h — * in (4.181 we get 



lim 

h^0 



Jn 



m [$~^ dxdt = lim j* ^ K6 [-L - I] ^ dxdt > 0, 



where the last inequality follows from Lemma 4.2 and the convexity of z 1— ► -. 
We thus have 



> a.e in SI x (0, T) , 



and the claim (4.151 follows from the lower bound [CI] on (j,(9). 
Step 2: By (14.141) and (|4T5| it follows that 



logT - logT = 



> 



a.e. in (0,T) x SI. 



(4.19) 



As log To — > log to a.e. in SI we conclude that 

logr (as, t) ~ log r(a;, t) > for a.e (x, t)efix (0, T) . 



On the other hand, Lemma 4.2 shows that \ogr{x,t) — logr(x,t) < a.e. Thus, 
logT = logT a.e. in SI x (0,T), and the conclusion follows by an application of 
the last part of Lemma |4.2| □ 
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4.3 Concluding the proof of Theorem |4.1 

In view of the two previous lemmas and [CI], 
\r h )^p(9,r) 



p(9 h ,r h )^p(9,r) inL 2 (0,T;L 2 (f})), 



rh 



in L p (!lx (0,T)), VI <p < oo . 



(4.20) 



Given <j> G C£°(fi x [0,T)) we use (3.5 1 with v h = DX [0] and integrate in time: 

cT 



u h Y\^ [<f>t] dxdt + / UqII^ y>o] dx 



F h n%[<t> x ] dxdt. (4.21) 



la Jn Jn J o Jn 

Due to the regularity of <f>, standard interpolation estimates show, together with 



[CO], (4.11, and (4.201, that the limit h -> in (4.211 results in (2.4 1. By now we 



have established the a.e. convergence (r , U , 9 h ) — > (r, u, 9), and that the latter 
triple satisfies (2.3| and (2.4 1. It remains to prove that (t,u,9) also satisfies 



(2.51. In order to pass to the limit in the temperature scheme, the following 

u, 9) s 

Then, 



lemma is essential; its proof exploits the fact that (t,u,9) satisfies (2.4 1 
Lemma 4.7 Let {r h ,u h ,9 h ) be as in Theorem 



4.1 



lim 

h^O 



sup 

te(o,T) Jn 



r r n(e n ) h a 

dx + / — ^tij-u, dam 
Jo Jn T 



0. 



Proof: Using w as test function in (3.5 1 and integrating in time we get 



l-of 



dx 



K0 h )^-p(0 h ,r h )uU^dt. (4.22) 



Sending ft, — > in (4.221, and using (4.201, gives 
i h (t)\ 



lim 

h^O 



+ 



M (6» h )i^| 2 ds r// 





114 


2 








77 




2 




Jo Jn 



On the other hand, by Lemma |4~4| u t G b L 2 (0, T; W _1,2 (f2)). This and the 
weak form of the momentum equation (2.4 1 yield the weak form 

rt 



[ ( u t, v ) (w -i 
Jo 



w 1,2 > 



valid for all v G 1/(0, T; W c , (fl)). Applying this with v — u gives 



IN')H! , r /•..^l«-l a <wt= W3 ' 



2 + .L.L^ 



+ 



whence 



lim 

h^O 



K*)II2 ll«o||l 



+ 



jj.(9) - — — dxdt 
r 



= Vt G (0, T) . 



Next, m£ ->■ u x in L 2 (0, T; L 2 (Q.)) and -» ^ in L P (Q,T; L p (tt)) for any 

p < oo. Consequently, in view of Lemma }4.1| 



lim f [ 

h ^oJo Jn 



^—r-^-u^.u x dxdt — 
Th Jo Jn t 



ui dxdt . 
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such that 



lira 

h->0 



sup 

t€(0,T) Jn 



- — — ^ dx+ f f ^ , \ux — u x \ 2 dxdt 
1 Jo Jn T 



= lim 

h^O 



sup 

te(o,T) Jn 



,/,■+ / T f /i{6 h )\-^f- - dxdt 



io Jn 

As the left hand side is non-negative this concludes the proof. 



□ 



Finally, fix any <j> £ C^°(Q x [0, T)). A calculation using (4.61 shows that 



i=l n "» 
;t, r,„ /Tq - ^ M™, rrQ 



((?")„ (n«M) B d!C, (4.23) 



where 0ft. is as in (3.91. Now using [(f)] as test function in (3.61 we obtain 

ri ft i 2 i 2i 

L-^J^drrdi. (4.24) 



3 h (f>t dxdt — I 0Q(l>odx- 



o Jn 

T 



\p(G h ,T h )u x 



On each interval (a;^ , x i+ x), n{6%) is a linear combination of k,{9^) and k(#^ +1 ). 
Hence, by Lemma |4~5| and since /3 < 2, «(#£.) -> k(9) in L 2 (0, T; L 2 (Q)). Thus 



4.1 



/^(a;))^)* ->■ k{9)9 x in i 1 (0,T; L 1 (f2)). Also, as ^ -> A in L p (0, T; L p (f2)) 
Vp 6 [1, oo), Lemma 



gives 



K(0) 



9 X mL 1 {Q,T;L 1 {Q)). 



Letting h — » in ( |4.24[ | and using Lemma |4~7| and ( |4.20[ ), results in ( |2.5[ ). This 
concludes the proof of Theorem |4.1| 



5 Energy bounds and proof of Theorem 2.1 



We now restrict to constitutive relations as described in Theorem |2.1[ for which 
the pointwise estimates ( |2.10[ | are established in the next section. Taking these 
for granted we now prove the integral bounds ( 2.11[ ). As demonstrated below 

an|2.1| 



these suffice, via Theorem |4.1[ to establish Theorem |2.1| For reference we note: 
[Al]: We have a = and (3 6 [0, |) such that T K and L are given by 



F = —i,u x -p{9 ,t 

T 



and L{9) = 



R9?+ 2 
( ( S + l)(/3 + 2) 



[A2]: {T h ,u h ,9 h ), h > 0, are solutions of the scheme (|3~i}-(|3~6t with initial 
data (tq ,ug ,6q ) given by (3.31 with (ro,uo,#o) as in Theorem 



2.1 



For a fixed T £ (0, oo) we let C, C, etc. be numbers that depend on T, system 
parameters (Q, fi, K), and the initial data, but that are independent of h. 

Lemma 5.1 (Pointwise estimates) Assume [Al] and [A2]. Then there ex- 
ists a number C > which is independent of h and such that: 



er 1 <r h (x,t) <c, v(i,t)efix(o,T), 



(5.1) 



14 



c 1 < e h (x, t), v(«, t) e n x (o, t), 



(5.2) 



lf h (*)ll. 



= («) 



dt<C. 



(5.3) 



The proof of ( 5.1 H 5.3 1 is essentially the same as in 16 . Minor adjustments 
are required to treat the particular scheme ( |3.4[ )-( |3~6| and to incorporate 9- 
dependence in the heat conductivity n. For completeness we include the proof 
in Section [6] We note that these pointwise bounds do not seem to generalize in 
any simple way to the case of ^-dependent viscosities. 



We proceed to state the discrete analogue of (2.11 1 (Lemma 5.2 below), and 
then show how this is used together with Lemma [5.1| to prove Theorem [27T] The 
remaining parts of this section detail the proof of Lemma |5.2| A few technical 
lemmas are collected at the end of the section. 

Lemma 5.2 (Energy estimates) Define 

A h (t) := sup / \iix\ 2 dx + / \ut\ 2 dxds, 
s£(o,t)Jn Jo Jn 

B h (t):= sup / L(9 h )dx+^J2 / G l {r h )lL'{9 h )f i ds, 
s6(o,t) Jn n ~[ Jo 



V h {t) 



i JV ~ 1 r* r 

su p ^rE Gi ( r ^ L '( e )]*+ / / 

se(o,t) ln ~[ Jo Jn 



K{6 h )\9t\ dxds. 



(5.4) 



Then there is a number C > 0, independent of h, such that 
A h {t) + B h (t)+V h (t)<C Vte(0,T). 

5.1 Proof of Theorem O 

We now take Lemmas |5.1| and |5.2| for granted, and we verify that these are 
sufficient to verify the conditions in Theorem 4.1 (with a = 0). First, our as- 



sumptions on the initial data in Theorem 1 2 . 1 1 are stronger than the corresponding 



con ditions in Theorem 4.1 Next, [CI] is an immediate consequence of Lemma 
5.1 [C2] f ollow s from the bound on B h together with the pointwise estimates 



of Lemma |5.l| and [C3] follows from the bound on A h . We can thus apply 
Theorem |4. 1| and conclude the existence and convergence parts of Theorem |2.1| 
At this point, Lemma 



5.1 



yields (2.10l 



It remains to verify (2.11|, and we consider the first term in the 2?(t)- 
functional in detail. By the T> h (t)-bound in (5.4 1, we have as in (4.231 that 



|tt(fft)| 2 



^frfi^E 6 '^)^)!'^ Vte(0,T)- (5.5) 



/ 

Jn 



In view of Lemma [4. 5| this estimate gives 

flfc.Aft, mL°°(0,T;L 2 (n)). 
Clearly, i, 9 6 L°°(0, T; L°°(f2)), such that 

sup f \L'(9) x \ 2 dx = sup f — - \9 X \ 2 dx 
te(o,T)Jn te(o,T)Jn r 



< C||0||i,oo( O T;I ,=o) — ll^llL^(0,T;L 2 (n)) < C. 

' II T ll£<»(0,T : £<»(n)) 



15 



The other terms in (5.4 1 are treated similarly and hence we conclude (2.11 1. □ 



5.2 Proof of Lemma 15.21 

Lemma 5.3 Assume [Al] and [A2] . Then there is a number C > 0, indepen- 
dent of h, such that 



sup 

te(o,T) Jn 



\u h (x, i)\ 2 dx + 



\u x (x, t)\ 2 dxdt < C. 



(5.6) 



Proof: Applying (13. 51 with v h = u h and fi = fi, yields 



(I 
dt 



f \^L dx+ f fl^a dx= f p (r h ,e h )u h x dx. (5.7) 
Jn z Jn T Jn 



Applying the Cauchy-Schwarz inequality with a suitable parameter e we have 



/ 

Jn 



p{T h ,d h )u h x dx < e [ Ar\v£\ 2 dx+ — 1| \ \\ L <*>(n) \\Q h \\t°°{n) 
Jo t e T 



dx . 



We choose e small enough that the first term can be absorbed on the left-hand 
side in (5.71. We then apply (3.81 together with (5.1 1. Integrating in time and 
using (5.31, yield (5.6 1 with a suitable C. □ 



Bound for A 



Applying (3.51 with the test function v h (x,t) := u^(x,t), and fj, = fi, gives 



|mJ"| 2 dx - 



flU x U tx h h , 

z P U tx dx 



2r h 



dx 



H\u x 



. h h j 

+ P u xt dx , 



where p h := p(6 h ,T h ). Integrating in time and using r t h = u x we obtain 



\u x 
2r h 



fT~(t) dx + / \ut\ dxds 



(5.8) 



2r h 



(0) dx 



pt — / h\3 r f ^ 

/ of h\2 +Pt u x dxds + / P h u x {s) dx\ =:V"ij, 
Jo Z \ T ) Jn i=1 



where we have also applied integration by parts to the pressure term. From the 
requirements on the initial data, we have that Ji < C. Next, by adding and 
subtracting the positive term \u x \ 2 p h /r h , we have 



Applying Lemma 



5.4 



o Jn 
with 



fj{n h x f 
(r h ) 2 



dxds < 



(u x ) T 



dxds 



fhy , gives 



-£-(x, s) dx 



< C |e [l + A h (t)] +lf* A h (s)\\u h x (s)\\ 2 LHn) ds 



1G 



Next, p h 



such that 



7 3 = ^/ /ML ^dxd S <C 



■A k (a)||0 h («)|Uoo (I1 ) ds + 2? h (*) 3 



where we have used (5.1 1 , (5.2|, and Lemma|5.3| To bound J4 we use the Cauchy- 



Schwarz inequality with a parameter and ( |5.2p together with the requirements 
on the initial data: 



I4 = I p h Ux(s) dx 
Jn 



< C 



sA h (t) + -B h (t) 



+ \\p h (0)\\L2(n)\\u x {0)\\ L 2 m . 



Using these bounds in (5.8|, taking the supremum over times in (0, t), applying 
( |5.1[ ), and choosing e sufficiently small, we obtain 



A h (t) <C\l + B h (t) + V h (t)i + J A h (t) [||t£| 



1,3(11) + II" I 



(it 



Applying Grdnwall's inequality together with (5.31 and Lemma 5.3 then gives 
A h (t) <c\l+B h (t)+V h (t)^. (5.9) 

Bound for B h 

To bound the B h functional we define the test functions 

^(M) :=L'(e h (x,t)), 

where L is given in [Al]. Using ip h as test function in ( |3.6| l, with fj, = fi, 
integrating in time, and rearranging yield 

r 1 JV_1 /•* 

' L(9 h (t)) dx+±J2 / Gi(r h )[L'(6l h )]]? ds 
a n j=1 Jo 

= f [ T h u x L\e h ) dxds + f L(8 h )(0) dx =: h + h ■ 
Jo Jn Jn 



Using the Cauchy-Schwarz inequality and the lower bounds in (5.1 1-(5.2 1, gives 

1 dxds + J* J \r h \ 2 \e h f +1 dxds] 
< c{ J* J \u h x f\e h \ B+ i dxds + J' J |0 h | a |0Y +1 dxds } 



h < a\j j \>i ! :Y\o 



< c 



{ [ jup \\e\s)\\f+* n) )] i + J* \\e h { s )\\ La o {n) [ J L(e h )dx] ds} 



C{l + V\t)i + J ||0»|U~ (n ) 



h (s) ds} , 



where Lemma |5.5| is applied in the last inequality. The I2 term is bounded by 
the requirements on the initial data, whence 



B 



h (t) < c[l + V h (t) 1 z + j \\0 h {s)\\L°°B h (s) ds} 



An application of GronwalPs inequality, where we use the bound (5.3 1, yields 

B h (t) < C[l +V h (t)?] . (5.10) 
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Bound for V h 

To bound T> h we define the test function ip h G Qji(f2) by 

i> h {x,t) = L'(9 h (x,t)) t = K(e h (x,t))9?(x,t). 



Using ip in the temperature scheme (3.6 1 and integrating in time give 



K{9 h )\e h t ? dxds+j R J2 J o G i (r h )(iL'(9 h )jl) t ds 
= f [ f h ui{L'(e h )) t 

Jo Jn 

Integrating by parts in time and rearranging gives 

= fJj h ul{L'{9 h )) t dxds + ^ E [Gi{T h )lL'{e h )t] (0) 
1 N ' 1 r* 



dxds . 



K,{8 h )\9t\ 2 dxds 



(5.11) 



Using the Cauchy-Schwarz inequality with parameter e together with the point- 



wise bounds (5.1l, (5.2 1 on r and 9 (and the assumption /3 > 0), we get 

Ji<e [ [ K(9 h )\9*\ 2 dxds + - I 

'o Jn e Jo Jn 

a r* , r r , . , . l 

dt 



K(9 h )\3 rh \ 2 \ux\ 2 dxds 



<eV h (t) + - f M9 h )\\ L o 
t Jo 



9 h \ 4 + \T h \ 4 dx 



l_L°°(Q) 



< eV h (t) + - \ sup / L(9 h ) dx 

H [ SU P IHOll^tn)! sup / \f h \ 2 dx] I \\T h \\ 2 L ^ (n) ds 

e L se(o,t) Ji se(o,t) Jn J VJo 



C . 



Jo 

Mo h )\\ L ^ m ]\. 

-se(o,t) 



ds 



G [ 



ds 



\ sup M9 h )\\ L oo (n) ]\A h (t) +B h (t) 
Applying Lemma |5.6| and the bounds ( |5.9[ ) and ( |5.10| l, we obtain 
Ji < eD h (t) + -B h (t) [l+6 h (i)i] 

+ ^ [ sup || /e(0 h ) || LO o {n) jA h (t) + B h (t)][l + .4 h (t) 2 ] 

< e2?*(*) + — [l + 23' 1 (*)ll[l+ sup ||K(^)|U- (n) l . 

6 L J L se(o,t) ' J 

To bound the last factor we apply Lemma |5.5| 

ii«(tf k )iUo. (n) = a B u P [ii^iif+y^^cfi+^^^l 



sup 

s£(0,t) 



18 



Thus, 



Jl < eV H (t) + 4(2/3+3) 



(5.12) 



10/3+9 



where we note that 4(2/3+3) 
requirements on the initial data 



< 1, since /3 < 



Next, consider J2: by the 



N-l JV-1 

Ja = i E < ^ [ E>o1?] < CUflo^H^co). (5.13) 

i=l i-1 

To bound J3 we first observe that for each 1 < i < JV, 

(G l (r h )) t = -iG,^) 2 [u£k +u*k +1 ] • 

Writing u£ = ^(.F ft +p h ) and using the positivity of # h together with both 
bounds in (5.1 1, we get 



(n) , 



where G > 0, 



and thus: 



1 JV_1 /•< 



J 3 = 



G 







£ G t (T h )lL'(8 h )j* 



i=l 
JV-1 



11^ 



c/.s 



< 



g \ sup, e ^wrf /vv 



-»e(o,t) j=1 

< CT^i)? [l+A h (t)iy B h (t)i < C \l + V h (t)i^ , (5.14) 
where we have used Lemma |5.6| together with the previously derived bounds on 



A n and B h . Substituting ( |5.12[ ), ( |5.13[ ), and ( |5.14[ ) into |5TlJ we obtain 

I)'' 
such that 



) h (t) < G [l+X^t) 5 ] 



with 5:= max { f, < 1, 



2T(t) < G. 



Recalling (5.9 1, (5.101 we conclude that there is a G > 0, depending on T, the 



systems parameters, and initial data, but independent of h, and such that 

A h {t) + B h (t) +V h (t) < G, for all t G [0,T]. 
This concludes the proof of Lemma |5.2| □ 



Technical lemmas used in the proof 

Lemma 5.4 Assume [A1]-[A2]. Then there is a C independent of h such that 



T h (x, s)4>(x, s) dxds 



(x, s)\ dx 



ds 



for all e > 0, t G (0, T) and <f> G L 2 {Q,T; L 1 (Q)) . 
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Proof: For each t € (0, T) define the test functions v h by 



v h (x,t)= / K<f]{y,t) 



(z, t) dz 



dy. 



Using v as test function in (3.5|, rearranging and integrating in time, give 

1 



T 6 dxds = 



mi 













/' 


1 <f> dx 




j J- h dx 


ds 


10 


Jn 




Jn 





n 

< eA h (t) + - 
e 



\(f>\ dx 



ds + C 

I 2 



dx 



(\u x \ + # h ) dx 



ds 



[/'< 


!>| dx 


Jn 





ds + e 



\Ux\ 2 dxds + Cet , 



where we have applied the Cauchy-Schwarz inequality with parameter e, together 
with the bound (3.8 1. Applying the estimate (5.61 concludes the proof. □ 



Lemma 5.5 Assume [A1]-[A2]. Then there is a C independent of h such that 



SUp H^f/o^ 



< C [l+V h (t) a } , 



sG(0,t) 

for all < q < 1. 

Proof: Fix any s £ (0, t). Since 6 h = 9 h (-, s) is piecewise constant, we have 

i-l 



9 h L'{9 h ) 2 \ 



L'{8 h ) 2 \ 



k=j 



= 6 h \ E L'(8 h f\ Ej ±0 h \ Ei Y,lL'{0 h )h[L'(0 h )\ Ek +L'{6 h )\ Ek J. 

Multiplying through by h, summing over j = 1, . . . , N, taking the maximum 
over i, and applying the Cauchy-Schwartz inequality, we deduce 



° h (s)l| 2 L l + ( 3 n) < ||fl h WIU- (n) f L'(9 h (s)) 2 dx 

j n 



2/3+3 



P\s)\\ 2 L ^ 2 n) + c\\e n (s)\\ L > in) 



N-l 



<49 h (s)\\f+^ )+ C(e) 



N-l 



where we have used the energy bound (3.8 1. We choose e suitably small, absorb 
the first term on the right-hand side into the left-hand side, and take the supre- 
mum over time. Recalling the definition of T> h (t), using the lower bound (5.1 li, 
and taking the ath power of both sides, then yields 



f sup \\e\s)\\f^ V <c\i + sup E \lL\e h ){s)fX <c\i + v h {tr] 

L se(o,i) 1 ' J L se(o,t) n J L J 



□ 
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Lemma 5.6 Assume [A1]-[A2]. Then there is a C independent of h such that 



f 

Jo 



\T h (s^" 2 



L°°(f2) 



ds<C\l + A 



and 



9 h (s)\\U m ds<C 



l (t)*] 



[l+B h {t)^ 



(5.15) 



(5.16) 



Proof: Consider the ith element Ei of the mesh and let denote the 
constant value of T h on this element. We have 

l.Fjf = i^f ±£ll^Tl*> 

k—j 

and we proceed to multiply by a^h, integrate in ti me, a nd sum over j. Making 
use of the lower bound on t' 1 , ( 3.8 1, (5.31 and Lemma 5.3 to bound J J \J- h \ 2 dxds 
we obtain 



/V* 

Jo 



| 2 ds<C 



\J rh \ 2 dxds 



< C{ 1 



< C< 1 



\J rh \ 2 dxds 



o Jn 

L 



h 2 



ds 



(5.17) 



where we have used the Cauchy-Schwarz inequality. To bound the remaining 
term, let Xj be an arbitrary interior node and let v G V®(Q) be such that 
v h (xk) = 1 and v (xi) = 0, VZ 7^ k. Using v as a test function in (3.51 we have 

IT T~hi\ I rr~h h 7 1 / T~h hi I h h 7 

W Ifc = - / Fk+\V X dy+ T k v x dy= u t v dx. 

J E k+ i J Ei. JBtUBtxi 



such that 



+ (n^Kl 



Since fc is arbitrary it follows that 



Y.iKf<c( {l%\u h t \Ydx<C [ |«?| 

fc — 1 



da; . 



Using this in (5.171 yields 



|^| 2 ds < C [l+A h (t)i] , 



and (5.151 follows. To establish (5.161 we argue similarly to get 

JV-l 



1\ 2 <\ [ \e h \ 2 dx+J2(e% + eZ +1 )\ie h ] 
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Using (3.81 and (5.31 to bound J J \9 h \ 2 dxds, together with (5.1 Ij, we obtain 



f\e1\ 2 ds<ch + J*J^(et + ei +1 )\ie h j h \ ds\ 

j E / G k (r h )lL'(e h )f k ds 



<C II 





i 




/ \9 h \ 2 dxds 


2 


[I. 


la, 







< C{ 1 + 



where we have used that \L' {8 h )\ k = K{8*)\8 h } k > C\8 h } k , (for a 8* between 
8 k and 8 k+ i) due to the lower bound (5.2 I. □ 



6 Proof of Lemma 5.1 



Throughout this section the assumptions [Al] and [A2] in Section [5] are in 
force, and we fix an arbitrary time T > 0. To simplify the notation we set 



h (x,t) := / T h (x,s) ds + IL% 



«o (y) dy 



In what follows C, d, . . . denote positive numbers that depend on T, \Q\, the 
initial data, and the bounds (2.9 1, but that are independent of h. 



Lemma 6.1 There is a number < Ci < 00 such that 

\io h {x,t)\ < Ci , for all (x,t) € X [0,T]. 
Proof: For fixed t > and v h S V®(Q.), integrating (|3.5[) in time gives 



h h 7 
LO dx 



h h 1 

u v dx . 



(6.1) 



Using v h (x) — u h (x,t) then yields 



Using this together with ( |3.4[ ) with 4> h (x) = cj h (x,t), give^] 
^ ./n 



h h j I rr~h h j 1 / h h j / r/- nh , I I 2 j 

w r dx = I J- t dx + I u> u„ ax = — I Jxa + \u dx . 



According to (3.8 1 we thus have, for each t > 0, that 



C -t(K + 2)8(9$, u%) < ( / TV(fe)(t)< / w' l (2;,0)r h (a;) dx =: C*. (6.2) 



Next, by conservation of mass ( |3.7[ ) there are maps i i— > z (i) g O such that 

^( 2 -(i),t) < M T-^ticJ (t) < uj h (z+(t),t), Vi. (6.3) 



x At this point the assumption of constant viscosity is used in an essential manner. 
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Next, fix t, h > 0, and y € Q arbitrary. Let E G Ed be the element containing 
z~(t), let E v be the element containing y, and define the interval S y C £1 by 

S a := [min{.z"(f),y},max{2~(t),j/}] U £T U E y . 

Then fix a v h g V^(n) by requiring that for each a^, i £ {0, 1, . . . , N}, 

1 iixi<ES v , 



v h ( Xi ) 



otherwise . 



Since io is piecewise constant we have 

u; h (y,t)-u; h (z-(t),t)=± l" h l< 

N-l 



± \^ [w Jiii (xi) = =F / w 8, dx = ± u h v h dx, 
i=1 in in 



where the last equality follows by (6.1|. Rearranging, applying Holder and the 



energy estimate (3.81, using ( 6.3 1 together with (6.21, and finally taking the 



supremum m y, give 



sup u) (y,t) < C . 
yen 



A similar argument using z + (t) yields inf H gn w (y, t) > — C. 



□ 



Lemma 6.2 There exists Ci > such that 

C2 1 < ||T h (t)||i«(n) < C* 2 [l+ / \\8 h (s)\\ L ~ {n) ds 

Jo 



for t € [0, T] . 



Proof: Since r/ 1 (■,£), «£(•,£) £ Q/^fi) at each time, ( |3.4[ ) gives r/ 1 = u£. Thus, 
with A := (r h ) _1 exp(/Z _1 Lj h ) and B := A' 1 , we get 



T, 



A < 0, 



S 4 



M M /" 



give 



Integration in time, together with (2.9 1 and Lemma 6.1 

A(x,t) < A {x) < C => r h {x,t) > C2 1 , 

and 



B(x,t) < B (x) + C I \\e h (s)\\ x ds T h (x,t)<C 2 



1+ / (s) || 00 ds 




□ 



Lemma 6.3 There exists C3 > such that 
1 



9*(t) 



<C 3 /orte[0,T]. 



(6.4) 



L°°(n) 
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Proof: Apply (3.6 1 with ip h 



-M(6 h 



(completing the square in u£) shows that 

N-l 



(i 
(It 



J^e h )- M dx-f ]T G^W) 



where M > 0. A calculation 



MK dx 



a 4/2^(^)^-1 



Since L' = k > and M > the sum on the left is negative, and Holder gives 



d_ 
dt 



1 

WW) 



MK 



dx 

r h(nh\M-l 



< CM 



where || • ||m denotes the L M (f2)-norm. Applying Lemma 



6.2 



1 

W) M 

we obtain 



d 
dt 



1 



< c 



T h (t) 



<c. 



Integrating in time, applying ( 2.9 1 2 , and sending M 1 oo, yield (6.4|. 



□ 



Lemma 6.4 There exists C± > such that 

rT 



ds<C 4 - 



(6.5) 



Proof: Fix i, j 6 !!. As 8 h (-,t) G Q^(f2), Lemma 6.3 yields (suppressing t) 
9 h (x)i = 0' l (y)e h (x)-5 + V < C^(y) + V ittl , 



where x € Ek and y £ Ej. Integrating in y, using (3.8 1, and choosing x such 
that \\e h (t)\\ L oo (n) = 8 h (x,t), yield 



'toll 



_L°°(f2) 



< c 



? h (*)ii 



L°°(n) 



< C 



where is as in (3.9 1. Multiplying and dividing by suitable terms in the latter 
sum, using k(9) oc 8 p , applying Cauchy-Schwarz, and squaring, give 



'toll 



< c{ 1 + 



E 



Gi(r")«(^ t ) Mfcl2 



h \2 



he?, 



E Gi{T h )(o^y 



(6.6) 



where f^j. is as in ( |3.9[ ). Let the two inner sums on the right-hand side of 
be 21(f) and ©(f), respectively. To bound 33(f) we recall that Gi{r h ) = 
2(r h | i3i + r h | Bi . +1 ) _1 , and use that 0& < 8 h \ E , + h \s i+1 , while 8^ is between 
8 h \Ei and S h |_E i+1 . Lemma 6.3 ( 3.8 1, and Lemma 6.2 then give 

»(*) < <5||r' i (t)|U ooCn) < c[l + jf ||fl h (a)|U- ( n) d«] . 
Squaring both sides in (6.61 thus gives 

z.<» ( n) <C 1 + 21(f)- [l 



'toll 



||0>)IU» (n) ds 



Applying the Gronwall inequality and recalling, by (3.10 1, that J Q 2l(s) ds < Co, 
we obtain (16. 51). □ 



Combing Lemma |6.1| - Lemma |6.4| completes the proof of Lemma |5.1| 
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